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Abstract 

The effects on the exchange coupling in magnetic trilayers due to the presence 
of a spin-independent potential well are investigated. It is shown that within 
the RKKY theory no bias nor extra periods of oscillation associated with 
the depth of the well are found, contrary to what has been claimed in recent 
works. The range of validity of the RKKY theory is also discussed. 
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Metallic magnetic multilayers are systems composed of alternating ultrathin layers of 
magnetic and non-magnetic metallic materials. They exhibit exchange coupling between 
magnetic layers across a non-magnetic spacer. Such a coupling oscillates between ferromag- 
netic and antiferromagnetic and its strenght decays as the spacer layer thickness is varied. 

Basically, two mechanisms have been proposed for explaining the oscillatory coupling in 
multilayers. Edwards et al. @|| have related the phenomenon to the existence of quantum 
wells for electrons (or holes) with both spin orientations propagating through the multi- 
layer structure, which result from the exchange interaction inside the magnetic layers. The 
existence of those quantum wells was lately confirmed by photoemission experiments ||. 
According to the model proposed by Edwards et al., the observed oscillation in the coupling 
results from quantum interference effects inside the well and bears a formal analogy to the 
de Haas- van Alphen effect. The other proposed mechanism for the coupling in multilayers 
is due to Bruno and Chappert and is based on an extension of the RKKY theory to the 
particular geometry of the system. In their approach the coupling arises from the polariza- 
tion of magnetic carriers in the spacer, as in the ordinary RKKY theory for the coupling 
between localized magnetic moments in a non-magnetic host. 

The formal connection between the two theories was recently made explicit by d'Albu- 
querque e Castro et al f|. In their work, model-independent closed form expressions are 
derived for the coupling, as well as for the bilinear and biquadratic terms, and it is shown that 
the full confinement and the RKKY-like theories correspond to limiting cases in which the 
exchange interaction inside the magnetic layers is very strong and very weak, respectively. 
Another important point made by d'Albuquerque e Castro at al. is that in the usual RKKY- 
like theories the difference between the non-magnetic parts of the potentials in the spacer and 
in the magnetic layers are neglected, which may lead to significant errors in the amplitude 
and phase of the coupling. 

Both the quantum well and RKKY-like theories relate the period of oscillation of the 
coupling to the extremal points of the spacer Fermi surface in the direction perpendicular 
to the layers. However, in a recent paper by Jones et al. || the question of other possible 
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sources for the period of oscillation is discussed and it is claimed that the oscillation period 
is determined by the depth of the quantum well formed by the magnetic and spacer layers 
rather than by the dimensions of the spacer Fermi surface. In addition those authors claim 
that such effects can generate a bias in the coupling, producing a prefered ferromagnetic 
coupling. In recent calculations Munoz and Perez-Diaz ]TT| have also found periods in the 



coupling not related to the spacer Fermi surface, which they claim are associated with the 
confined states in the quantum well. 

In this communication we present results for the coupling for cases in which the quantum 
well due to the difference in the potentials of the spacer and the magnetic layers give rise to 
bound states for carriers with both spin directions. Those are just the situations considered 
by Jones et al. and Munoz and Perez-Diaz. We have found that, provided the contributions 
to the coupling coming from singularities associated with the presence of the bound states 
are properly dealt with, no extra period nor any bias are obtained. We therefore conclude 
that the features observed by Jones et al. and by Munoz and Perez-Diaz are spurious. 

We consider trilayer systems described by the one-band tight-binding model with nearest- 
neighbour hopping t on a simple cubic lattice. The layers are displaced perpendicularly to 
the (100) direction and the site energies are chosen equal to 6\t\ in the spacer and 6|i| + V 
in the magnetic material. Such a choice of parameters places the bottom of the spacer band 
at the origin of the energies. For sufficiently small values of the potential barrier V and 
of the Fermi energy Ep, the Fermi surfaces of both the spacer and magnetic materials are 
nearly spherical. Therefore the effective mass approximation can be safely applied to both 
materials, and the present model becomes equivalent to the electron gas model. We follow 
Jones et al. and introduce a local exchange interaction V ex only in the magnetic atomic 
layers next to the spacer. We label these two layers and n, so that the number of atomic 
layers in the spacer is equal to N — n — 1. The restriction of the local exchange interaction 
to just one atomic layer in the magnetic material would have to be lifted if we were to 
study the dependence of the coupling on the magnetic layer thickness. A recent thorough 
investigation of such a dependence has shown that as far as the behaviour of the coupling 



as a function of the spacer thickness is concerned, changes in the magnetic layer thickness 
just leads to changes in the phase and amplitude of the coupling, keeping the period of 
oscillation unchanged. Therefore, for the purpose of the present work the restriction of V ex 
to just one monolayer in the magnetic material does not pose any limitation on the validity 
of our final conclusions. 

The change in the thermodynamical potential Q at T = due to a rotation by an angle 
9 of the magnetization in one of the magnetic layers relative to that in the other layer is 
given by f§ 

An(6) = [ Ef dLuImln{l - 2V£ ( cos 9 - 1 ) G^u, gj,) Gj>, q {] ) } (1) 

7T _ J — oo 

where G^(u,qu) is the off-diagonal matrix element between sites and n of the retarded 
Green's function for an electron with spin a in the ferromagnetic 9 = configuration of the 
system. Here the summation over gj| is restricted to the two-dimensional Brillouin zone. It is 
worth stressing that the Green's function G a properly takes into account the presence of the 
potential well and, consequently, exhibits poles at the corresponding bound states energies. 
Those poles lie on the real axis. However, because the integrand in eq.(|l]) is analytic in the 
upper-half of the complex energy plane, the integration over the energy uj can be performed 
along a straight line joining the points Ep + ioo and Ep, thereby avoiding the singularities. 

The expansion of the integrand in eq.(|l|) in powers of cos(9) || enables us to introduce 
the bilinear J\ and intrinsic biquadratic J2 coefficients, in terms of which the expression for 
AO reads, 

AQ(0) = Afi - Jicos9 - J 2 cos 2 9 (2) 

It is found that higher order terms than cos 2 9 are negligible so that 

Ji = Afi(vr)/2 (3) 

Here we are interested in those situations in which the local exchange interaction V ex is 
small and the coupling can be calculated by perturbation theory in that parameter. It is 
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then straightforward to show that to lowest order in V ex the bilinear coefficient is given by 



J? KK¥ = - E f'dulm {2V^G° ) (u,q\\)GS n (u,q^)}, (4) 

7T -. J — oo 
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where the Green's functions are now calculated for V ex = 0. This is just the second-order 
perturbation theory or RKKY result for the coupling. 

Using the above expression we evaluated J^ K for the trilayer system as a function 
of the spacer thickness N, for fixed values of Ep and V ex , but for different well barriers V. 
For the systems under consideration it is possible to derive analytical expressions for the 
off-diagonal Green's function elements G^ (w, q\\) and GQ n (uj, q\\) for an arbitrary value of n 



T4 . Those expressions can be analytically continued to non-integer n, which permits us 
to consider continuous variations of the spacer thickness. We set Ep = \t\, which gives rise 
to a period of oscillation of exactly three atomic layers, and V ex = 0.04|t|. In all cases, the 
calculated values of J± KKY (diamonds) are compared, following eq.(|), with Af2(7r)/2 (solid 
circles), the latter quantity being calculated from the full expression ([!]). 

Fig. ([TJa) shows results for the case in which V = (no potential well). As expected the 
coupling oscillates about zero (no bias) and with a well defined period of three atomic layers. 
However, in contrast with what was found by Jones et al. and Munoz and Perez-Diaz, those 
features are not changed when non-zero values of V are considered. We have performed 
calculations for V equal to 0.3Ep, 0.6Ep, and 0.9Ep- Results are shown in Figs. (|l]b), 
(|I|c), and (|l]d), respectively. The only noticable changes in the curves are in the phase and 
in the amplitude of oscillations. We have observed, however, that because of the singular 
behaviour of the integrand along the real energy axis, spurious features, such as bias and 
extra period, can be easily obtained when the contributions from the singularities are not 
properly taken into account. We have explicitly verified this point by evaluating the energy 
integration also along the real axis. It should be stressed that numerical problems may occur 
even when V = 0. In such a case, for each the integrand in Eq.(4) has a pole right at 
the bottom of the band. In all cases, if the contribution from poles is not treated properly, 



a bias in the coupling is obtained. This point was already noticed by Yafet jn| in the case 
of the RKKY interaction for the one-dimensional electron gas with V = 

It is interesting to notice that for the value of V ex considered above, the second order 
perturbation result J^ KKY for the coupling agrees well with the full calculation AQ(ir)/2. 
For higher values of this parameter, however, the second order perturbation theory does not 
provide a satisfactory description of the coupling and full calculations have to be performed. 
This point is illustrated in Fig. (||D which shows results of both J^ KKY (diamonds) and full 
calculation (solid circles) as functions of V ex , for different values of V and N. We clearly see 
that for fixed Ep the range of applicability of the second order perturbation theory strongly 
depend on both N and V . 

In summary, we have shown that the inclusion of interface effects through a spin inde- 
pendent square well potential does not alter the period of oscillation of the RKKY bilinear 
coupling, as it was recently claimed. We have also checked that the appearance of a ferro- 
magnetic bias is due to numerical inaccuracy in the calculation and does not reflect any real 
effect due to presence of bound states. Finally, we have explicitly shown the limited validity 
of the RKKY-like theories. 
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FIGURES 

FIG. 1. Exchange Coupling as a function of spacer thickness for barrier heights V = (a), 
0.3Ef (b), 0.6Ef (c), and 0.9Ef (d). In each graph, results for J^ KKY are represented by diamonds 
and those for AQ(tt)/2 by solid circles. Full lines are for continous values of n. 

FIG. 2. Exchange Coupling as a function of the local exchange interaction V ex for different 
heights of the well and spacer thicknesses. Results for J^ KKY are represented by diamonds and 
those for AQ(tt)/2 by solid circles. 
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